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The generation of dark solitons in Bose-Einstein condensates with phase imprinting is studied by 
mapping it into the classic problem of a damped driven pendulum. We provide simple but powerful 
schemes of designing the phase imprint for various desired outcomes. We derive a formula for the 
number of dark solitons generated by a given phase step, and also obtain results which explain 
experimental observations. 
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PACS numbers: 05.45.Yv, 03.75.Fi, 42.65.Tg 

Solitons have been discovered in various classical non- 
linear media, such as fluid, magnetic and optical sys- 
tems, and have fascinated physicists for decades for their 
particle- like properties Q. Recently, dark solitons were 
observed in Bose-Einstein condensates (BECs) of dilute 
atomic gases ^ |!| , which are described by a macroscopic 
wave function. Dark solitons are produced by engineering 
the phase of this wave function with a technique known 
as phase imprinting, which was originally proposed and 
used to create vortices Q. Phase imprinting is to shine 
an off-resonance laser on a BEC thus create phase steps 
between different parts of the BEC. As a new tool of ma- 
nipulating a wave function, its power and ability have 
not been studied in a systematic manner. 

In this Letter we present a thorough analysis of 
the generation of dark solitons with phase imprinting 
in BECs. Our study is facilitated by a novel ap- 
proach, which maps the soliton generation problem into a 
damped driven pendulum problem. This method makes 
it easier to find the dark solitons generated by a given 
phase step. More importantly, it changes our perspective 
on the problem of soliton generation. With this method 
we can show how to design and control the phase steps for 
various desired outcomes, such as a specific dark soliton 
and a specified number of solitons. We derive a formula, 
which relates the winding number of the pendulum mo- 
tion to the number of dark solitons generated by a given 
phase step. In addition, we study the interesting phe- 
nomenon that counter-propagating dark solitons may be 
generated by one phase step, as observed in Ref.Q, and 
the physics behind it. Although our study is done in the 
context of a BEC, it can easily be applied to fiber op- 
tics, where dark solitons have potential applications in 
communication ||. 

We study the phase imprinting on a quasi-one dimen- 
sional BEC, which is realizable experimentally p] and 
was indeed used to produce dark solitons in BECs ||. 
On the other hand, although the magnetic trap is im- 
portant for the dynamics of BEC dark solitons R] , it has 
negligible effect on the generation of dark solitons, which 
has a much shorter time scale than the subsequent dy- 



namics. Therefore, it is sufficient to neglect the trap and 
use the one dimensional nonlinear Schrodinger equation 



.dtp(x, t) 



1 d 2 ip(x,t) 
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ul\iP{x,t)\ 2 iP(x,t) , (1) 



where x is measured in units of £ = 1 /im, a typical length 
unit in this type of experiments, t in units of — (m is 
the atomic mass), ip i n units of the square root of no, the 
average density the condensate, and the speed of sound 
is Mo = v 47rnoa s £ 2 , with a s > being the interatomic 
scattering length. For the experiment with rubidium ||, 
we have uq ~ 5.5 and the time of evolution around 10 ; 
for sodium experiment we have uq ~ 1 and the time 
of evolution around 30 . 

A dark soliton is characterized by a local density mini- 
mum moving with constant speed against a uniform back- 
ground H ||. It has three characteristics, the depth of its 
density minimum, the phase step over its density notch, 
and its velocity. However, all the three are related to each 
other and can be specified by its velocity. The nonlinear 
Schrodinger equation is exactly solvable with the inverse 
scattering method [|| , according to which the generation 
of dark solitons is determined by the Zakharov-Shabat 
(ZS) eigenvalue equations, 



. dUi{x) 
dx 

. 8U 2 (x) 
dx 



u ip(x,0)U2(x) = \Ui(x), 



(2) 



uo1> m (x,0)Ui(x) = -XU 2 (x), (3) 



where ^(x,0) is the initial condition. The ZS equations 
can have discrete eigenvalues with magnitude smaller 
than uq. Corresponding to each Aj, a dark soliton with 
velocity —Xi is generated. For phase imprinting, we have 
tp(x,0) = e lS ( x \ where S(x) is the imprinted phase. In 
this letter, for simplicity, we will concentrate on the right 
phase step, which increases monotonically from the left 
to the right, and stays constant at the boundaries, 



dS(x) 
da; 
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dS(x) 
dx 
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It is straightforward to generalize our method and results 
to the left phase step and more general phase imprinting, 
or even to the density engineering || . 

We solve the ZS equations by mapping them into a sim- 
ple pendulum problem, which is physically more intuitive 
and mathematically much simpler. Note that, for dis- 
crete eigenvalues Ai, the ZS wave functions U\ = ^fple 1 ^ 1 
and U2 = y/f>2e 1 ^ 2 have the boundary conditions 

/°i 5 /°2^0, 4>i , 4>2 — > constants , as |x| — > 00 . (5) 

As is well known, the quantity \Ui\ 2 — IC/2I 2 is conserved 
and independent of x, which, combined with Eq.(j^), leads 
to the conclusion that |J7i| 2 = \U2\ 2 = p for a discrete 
eigenvalue A,-. In light of this, we make the following 
transformation 

Ui = i^-pe^-^' 2 , U 2 = Vpe^- S )/ 2 , (6) 

which turns the ZS equations into a pair of very simple 
equations 

tp= ( -^- = 2\ + S -2u a smtp, (7) 
da; 

P = -7- = 2uo/ocos(^, (8) 
da; 

with 8 = 0, where the overhead dot denotes the spa- 
tial derivative. Remarkably, equation (0) involves only 
ip, and can be viewed as a damped massless pendulum 
driven by the force 2A + S if we regard x as time. This 
pendulum has two fixed points, P s and P u , when S = 
and |A| < uq. The point P s is at ip = sin -1 ^- and is 
stable; the other one P u is at 7r — tpo and is unstable, as 
shown in Fig . 

The solutions of Eq.(^) transformed from the ZS eigen- 
functions for a discrete eigenvalue always start at P s and 
end at P u . This can be checked by noticing that the 
boundary conditions (|J) become p — > and ip — > 0, and 
the discrete eigenvalues have magnitude smaller than uq. 



The correspondence between the pendulum solutions 
going from P s to P u and dark solitons can be appreci- 
ated in context of the pendulum problem (Q) itself with- 
out referencing to the ZS equations. With a given phase 
step, the pendulum equation (Q) has a solution starting 
at the stable fixed point P s for each |A| < uq, as shown 
in Fig .|l|. For future convenience, we name this kind of 
pendulum solutions proto-soliton solutions. Because of 
the asymptotic behavior, S = at |x| — > 00, the proto- 
soliton solution can only end up at either P s or P u . Since 
P s is the stable fixed point, for most values of A the pen- 
dulum comes back to P s after several rounds of rotation 
(Fig .0A). Only for finite number of A's the pendulum will 
end up at the unstable fixed point P u (Fig.^B). When 
this happens, we say the proto-soliton solutions become 
the soliton solutions. This set of special A's are just the 
discrete eigenvalues Ai of the ZS equations while the soli- 
ton solutions correspond to the ZS eigenfunctions. 




-71/2 71/2 371/2 27C 571/2 371 77C/2 47t 9tc/2 

9 



FIG. 2: Phase trajectories of a damped driven pendulum of 
zero mass. Dotted line are drawn for constant forces 2A; solid 
lines are for changing forces 2A + S. The arrows indicate the 
directions of pendulum motion. The black dots are the stable 
fixed points P s ; the open circles are the unstable fixed point 
P u . 




FIG. 1: Motions of Pendulum ([7|). Trajectories are schematic, 
and deliberately distorted when the overlapping occurs. The 
vectors are forces. A) Motion starting from P s and coming 
back to P s after one round of rotation. B) Motion going from 
the stable fixed point P s to the unstable fixed point P u . 



This novel approach has tremendous advantages over 
the existing methods for the study of ZS equations flio|| . 
With the above analysis it is clear that we can dis- 
card Eq.(||) and focus only on the pendulum equation 
(0), which is much simpler than the ZS equations. As 
an example, we solve Eq.rtTl) for the case of the "sud- 
den" limit in which the phase imprinted is a step func- 
tion, S(x) = a<d(x). In this case, the force is a 8- 
function, S = a8(x). Integrating Eq.(^), we have 
¥>(0+) — <p{0-) = (it — ipo) — (fa = a which gives us 
7r — 2<po = a. So only one soliton is generated, and its 
speed is A = uosin^o = uoCOs(a/2). This recovers the 
result in Ref.|I]. 

More importantly, this method changes our perspec- 
tive on the problem of generating dark solitons. It allows 
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us ask and answer the inverse question, "what phase step 
is needed to produce a specified dark soliton?" This is 
achieved with the following steps: 

1) pick A, the soliton that one wants to create; 

2) choose a curve, ip — f(ip), connecting the pair of fixed 
points P s and P u corresponding to A in Fig .0 

3) substitute ip — f(ip) into Eq.(^), solve it for S(x). 
The obtained phase step S(x) creates the soliton A. Note 
that any curve lying entirely in the upper half of the plane 
correspond to right phase steps, S > 0. 

Obviously, there are infinite number of paths connect- 
ing a pair of P s and P u , thus there are infinite number of 
phase steps that generate a certain dark soliton. We seek 
one best step in terms of expenditure of energy. Imprint- 
ing a phase step S(x) injects energy into the system, 
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4A - 4w sin tp + 2S 



(9) 



The minimum value of this energy corresponds to the 
smallest disturbance to the system by phase imprinting or 
also the smallest amount of noise in the output. Viewing 
E as a functional of S, using variational analysis, we have 
the phase step of the least disturbance, 



S(x) = 4 tan 



u - A 



tanhOrJug - A 2 ) . (10) 



This phase step generates a soliton of velocity —A, and 
possibly other solitons. Curve 1 in Fig .^| is one example 
of this kind. 

One interesting aspect of the phase step (|l^) is that the 
step height increases as A gets smaller when A is positive. 
This agrees with the experimental results that the soli- 
ton speed decreases when the step height increases [^|, || . 
More interesting is when A becomes negative. For A < 0, 
our numerical study shows that the step ( |I(i| ) always gen- 
erates more than one solitons with one traveling to the 
right with speed |A| and the rest to the left. That im- 
plies that counter-propagating solitons can be generated 
by one phase step, as observed experimentally ||. This 
phenomenon is certainly not special to the phase step 
(|l0|), it happens to any phase step obtained from a path 
connecting a pair of the fixed points of a negative A in the 
upper half plane of Fig However, this phenomenon is 
rather mysterious and against physical intuition. 

Intuitively, once a right phase step (ji|) is imprinted on a 
BEC cloud, atoms in the imprinted areas will start mov- 
ing to the right. Since the atoms outside of the imprinted 
area do not move, a dip with a bump to its right will ap- 
pear as a result (see Fig ||) . Due to the stronger repulsive 
interaction from the bump, the dip will be pushed to the 
left. As dark solitons come into form in the dip, one 
would expect that the dark solitons generated by a right 
phase step would always move to the left. Our numer- 
ical simulation shows that this intuitive picture is only 
correct at the beginning of time evolution. But after a 
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FIG. 3: Generation of counter-propagating dark solitons. 
Two snapshots are taken at t = 3.0 and t = 4.0. This initial 
phase step is a tanh function, S(x) = a(tanh(2x/a) + l)/2 , 
with uq ~ 1.0, a = 3.6, and a = 3.57T. 



very short time, as shown in Fig one soliton in its 
early forming stage may change its direction and start 
moving to the right. More interestingly, the phase dif- 
ference AS" = S(+oo) — S(— oo) is reduced by 2ir when 
this happens, as seen in Fig .|3|C. This means that the ini- 
tially imprinted the phase difference AS* can not always 
be maintained in the course of time evolution [[l2| . 

We are also able to design a phase step which gener- 
ates a specified number of solitons. This is based upon 
the formula for the number of solitons generated by a 
given phase step. To derive the formula, we need to 
examine the winding number I'F(A) of the proto-soliton 
solution, which is defined as the number of rounds of ro- 
tation that the pendulum makes in the solution. For a 
given phase step S(x), we have W(X) < W(X') if A < A', 
and W(X') — VF(A) = 1 when there is only one eigen- 
value Xi between A and A' |l4|]. Therefore, the number of 
soliton generated by S(x) is N s = W(uo) — W(— u a ), 
and the number of solitons traveling to the right is 
N r = W(0) — W{— uq) . As an example we consider a 
special but still quite general case, a phase step which 
can be cut into pieces so that S can be considered as 
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constant within each piece. Finding the winding number 
for each piece and adding them together, we have 



"• cs,NT (£vV? + "^) 

- int (/ + v\/t-«» s ) + 1 - (11 > 

where the second integral is done over the intervals on 
which S > 4uo, and INT is the integer function which 
keeps only the integer part of a real number. 

To imprint a phase which generates exactly n solitons, 
we only need to find a phase step that has W(uq) = n 
and W(— u ) = 0. This is achieved by drawing a path 
connecting tp = n/2 and (p — 2nn + n/2, the fixed point 
of the pendulum driven by a constant force 2uq. At the 
same time we make sure it lies under the curve, <p = 6uq — 
2uq sin tp, the darkened dashed line in Fig The phase 
step obtained from this path certainly has W(uq) = n. 
Also this phase step has W(— uq) = 0. To see this, note 
that the phase step satisfies < S < iuQ. Therefore, we 
always have <p < at cp = n/2 for the case A = — uq thus 
the pendulum can never pass the position n/2 and make 
a full round of rotation. The simplest example is curve 3 
in Fig J5[ which represent a linear phase step generating 
exactly two solitons. 
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FIG. 4: A. Number of Solitons. B. Number of solitons 
traveling to the right. 

Finally, we apply our results to a simple but very useful 
case, the linear phase step: S = a/a for I as j < a; for 
|ar| > a/2. In the present experiments ffl, one adjusts 
the step height and the step width to create different 
phase steps, which can be well modelled by the linear 



phase step. For the linear step, Eq.([ll|) is exact, and it 
is used to compute the number of solitons generated for 
different step heights a and different step widths a, as 
shown in Fig .|]A. Also for the simple case, Eq.(Q) can be 
solved analytically, and the eigenvalues Ai are given by 

A 2 -^ + tr _ >fo + fe) 2 -"3 (12) 

V u o - A2 tan(a</(A + - u 2 ) ' 

By solving this equation, we obtain the soliton velocities, 
and find the number of solitons traveling to the right, 
which is plotted in Fig .|JB. As a whole, Fig .[| serves as a 
reference table, where experimentalists can find the right 
parameters, the step height and step width, to generate 
desired number of solitons. 
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